NORDITA-2012-57 
UUITP-20/12 



Large AT Limit of Af = 2 SU(N) Gauge 
Theories from Localization 

J.G. RussoMElaad K. Zarembo 3 >f] 

1 Perimeter Institute for Theoretical Physics, 
Waterloo, Ontario, N2L 2Y5, Canada 
2 Institucio Catalana de Recerca i Estudis Avangats (ICREA), 

Pg. Lluis Company s, 23, 08010 Barcelona, Spain 
3 Nordita, Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden 
^Department of Physics and Astronomy, Uppsala University 
SE-751 08 Uppsala, Sweden 

j orge . russo@icrea . cat , zarembo@nordita . org 



Abstract 

We study J\f = 2 Yang-Mills theory on S* 4 in the large- N limit. We find 
that on a large sphere Wilson loops obey a perimeter law and that the 
free energy grows quadratically with the radius of the sphere. We also 
comment on the large-iV limit of the A/ = 2* theory, and on the free 
energy in J\f = 2 and M = 4 superconformal theories. 
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1 Introduction 

In supersymmetric field theories a certain class of non-trivial dynamical quantities can be 
computed exactly. An example, that we will further explore here, is the partition function of 
an Af = 2 supersymmetric gauge theory on S 4 . By using localization, the S* 4 partition function 
can be reduced to a finite-dimensional matrix integral pQ. Observables, that can be computed 
by this method, include the free energy and expectation values of circular Wilson loops. 

One can envision a number of applications of this remarkable result. The 't Hooft's large- N 
limit of the matrix integral [2] [3] H] (5J [6J may give us a glimpse on the putative string duals 
of Af = 2 theories. In particular, for the Af = 4 super- Yang-Mills theory (SYM), localization 
confirms earlier conjectures on the circular Wilson loop [7], which were instrumental in verifying 
the AdS/CFT correspondence. Or one can address generic questions about quantum field 
theory, such as the asymptotic behavior of perturbative series [6J in the setting where exact 
non-perturbative results are available. 

We will study the large-Y limit of the Pestun's matrix model, focusing on the simplest 
example of asymptotically free gauge theory, pure Af = 2 SYM. We will also comment on 
the Af = 2* theory (a massive deformation of Af = 4), which was studied in [BJ, and on the 
superconformal theories considered in [2] [3l HJ [5] . The motivation here is the AdS/CFT duality, 
and we will fill the gap in the literature by computing the finite part of the holographic free 
energy in Af = 4 SYM. 
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What do we expect to find, when we study an asymptotically free QFT on S" 4 ? In absence 
of other parameters in the problem, the dynamics is controlled by the running coupling renor- 
malized at the inverse size of the sphere (-R -1 ). When the radius is small, perturbation theory 
should work and yield an expansion of observables in inverse powers of lni? -1 . For a large 
sphere, the theory enters a non-perturbative regime, and the behavior of observables becomes 
a non-trivial dynamical question. We know that M = 2 SYM is not confining [S], and may 
thus expect that Wilson loops will obey a perimeter law, rather than an area law. As for the 
free energy, it cannot be extensive (cannot scale as R 4 ), because the vacuum energy is strictly 
zero in any supersymmetric theory. Since the perturbative free energy quadratically diverges 
in J\f = 2 SYM (see footnote 1), we may expect to find a A 2 R 2 behavior, where the UV cutoff 
is replaced by the dynamically generated mass scale A. The analysis of the matrix model will 
confirm these expectations. 

2 AT = 2 SYM 

2.1 Partition function 

The partition function of any M = 2 supersymmertic gauge theory on S" 4 is calculable by 
localization methods pp. The end result of this calculation is a finite-dimensional integral over 
the moduli space of Coulomb vacua, parameterized by the eigenvalues of the scalar field $ from 
the vector multiplet: $ = diag(ai, . . . a N ), tr $ = 0. We will mostly concentrate on the simplest 
example, pure M = 2 Yang-Mills theory without matter multiplets. Its partition function is 
given by the eigenvalue integral 

/ 8ir 2 N y a 2 „ 
d^aYli^-^fHH^-a^e A * 1 \Z- mst (a; g 2 )\ ■ (2.1) 
i<j 

The key point in the localization method is the one-loop exactness of the effective action for 
the eigenvalues. It thus suffices to compute one-loop determinants and to sum over instantons. 

The one-loop factor was calculated in pQ and is expressed in terms of the function H(x), 
defined as an infinite product over spherical harmonics: 

tfOO^n 1 + e~. (2.2) 

n=l ^ n I 

The coupling constant in the partition function is actually the running coupling renormalized 
at the scale set by the radius R of the sphere, the only parameter in the problem. Of course, the 
original path integral before localization depends on the bare coupling Aq, but this combines 
with the one-loop UV logarithms to produce the renormalized running coupling A. Indeed, 
the product over spherical harmonics appears in the calculation of [1] without the suppression 
factor exp(-x 2 /n) and thus logarithmically diverges. The infinite contribution, cut off at some 
n = A R is absorbed by renormalizing the bare coupling: 



as can be easily seen from (2.1) (7 is the Euler constant). A more mathematically precise way 
to get rid of divergences, used in [1] and also discussed later in section |4j is to consider M = 2* 
theory as a UV completion of pure M = 2 SYM. The mass of the adjoint hypermultiplet then 
plays the role of the UV cutoff A . 

The running coupling can be expressed through the radius of the sphere and the dynamically 
generated finite scale A = A exp(-47r 2 /A + 7): 

4tt 2 

= -In (AH). (2.4) 

A 



The running coupling becomes negative at R > A 1 . What makes the integral in (2.1 ) convergent 
is the rapid decrease at infinity of the one-loop function H(x). 

Although it is commonly assumed that instantons are unimportant at large-A, this is guar- 
anteed to be true only at weak coupling. In principle the entropy (the volume of the instanton 
moduli space) may overcome the exponential suppression by the instanton action. Then the 
instanton weight becomes exponentially large, and the system undergoes a transition to a new 
phase. Whether such transition happens or not is a dynamical question [9]. Such an instanton- 
induced phase transition does not occur in the Af = 2 superconformal theory [3J. This dynamical 
question is more intricate in the M = 2 SYM case, where the coupling runs and it is not clear a 
priori which is the relevant scale for A appearing in the instanton weight. As explained below, 
the scale at which the coupling is renormalized, can get dynamically shifted from R^ 1 -> aR' 1 . 
The new scale effectively represents the inverse of the instanton size, and a relevant question 
is whether the partition function is dominated by small or large instantons. These issues will 
be addressed in section 5, where the instanton factor |^i ns t| 2 in the large N limit is discussed. 
Meanwhile, we will first set 2i n st to 1, and in section 5 we will argue that indeed it is not 
relevant for the large N dynamics by explicitly computing the one-instanton contribution. 

Supersymmetric observables map, upon localization, onto expectation values in the matrix 
model. For example, the Wilson loop in the fundamental representation on the big circle of S* 4 
maps to the trace of the matrix exponential: 



W(C chcle ) = {fe 2 ™*). (2.5) 



Our goal is to solve the matrix model in the large-A limit and to compute the free energy and 
the Wilson loop vacuum expectation value (VEV). 

2.2 Toy model: SU(2) 

Before proceeding with the large- A calculation, it is instructive to consider the SU(2) partition 
function, which is given by a simple one-dimensional integral: 



yk=0 
Zj SU(2) 



-too 

J daa 2 H 2 (2a)(AR) 8a2 , (2.6) 
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which we will compute in the saddle-point approximation. As we will explain, the saddle-point 
approximation is not arbitrarily accurate for SU(2), but numerically works pretty well. A 
more serious problem is that one cannot really ignore instantons at small N, so this is just a 
toy model that illustrates the effect of the one-loop determinant. It gives a correct qualitative 
picture of what happens in the SU(N) theory at large N. 



The saddle point of the integral (2.6) is determined by a minimum of the effective potential: 

S c s(a) = -8a 2 \n(AR)-2\nH(2a) . (2.7) 

Here we ignore a term - In a 2 coming from the Vandermonde determinant, which can be viewed 
as a quantum prefactor. It is not significant for the large AR behavior. 

The function H(x) has the following large and small x asymptotics: 

, rr/ N f-^X 4 + 0(x 6 ) atS^O 

^h(x) = \ 2 y ; (2 8) 

-arm|:r|e' 2+0(ma;) at x ->■ oo 

The coupling constant here is renormalized at the inverse radius of the sphere R _1 , but only 
at small a. When a is large, 21niJ(2a) behaves as -8a 2 In a and a new scale aR^ 1 emerges. 
The effective potential thus has the Coleman-Weinberg shape [10], illustrated in fig. [I] The 
potential has a minimum at zero for RA < 1 (which is slightly shifted off zero if the Vandermonde 
determinant is included). However, if R > A -1 , this point becomes a maximum, and the saddle 
point shifts to a solution of the equation: 

AaJnAR = K(2a,), (2.9) 
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where 



K( x ) = -EQ = 2x £ f- - -3^-3 ) = x (V (1 + ix) + 1/> (1 - ix) - 2^(1)) , (2.10) 

jj n=l ^ ^ Tl ~t~ X / 

and ^(l) = _ 7- The saddle-point equation is essentially the condition g^a^Rr 1 ) - 00. For large 
radius, the saddle-point position a* thus grows approximately linearly with R, 2a* « e~ 7 Ai? 
(the saddle-point approximation applies when the second derivative of the potential is large; in 
the present case, S'J s (a*) = 16 + 0(R~ 2 )). 

The VEV of the circular Wilson in the saddle-point approximation is given by W(C C i rc \ c ) = 
cosh 2na*. Since a* = at R < A -1 , the Wilson loop stays constant at small R, within the 
semiclassical approximation. In practice, the exact VEV slowly grows with R until R reaches 
R c ~ A" 1 , after which a t *0 and the Wilson loop starts growing very fast. Asymptotically, at 
very large R, 

lnW(C c irde) xne-VAR, F = S(a*) - -e~ 27 A 2 i? 2 . (2.11) 

The Wilson loop thus satisfies the perimeter law, and the free energy is quadratic in the radius 
of the sphere. We will find the same behavior in the SU(N) SYM at large N. 

The perimeter law means that the mass of a heavy probe, associated with the Wilson loop, 
acquires a finite negative renormalization. The quadratic dependence of the free energy on 
the radius can be interpreted as a finite renormalization of the Newton's constant. In fact, a 
short-distance, quadratically divergent contribution to the Newton's constant does not vanish 
in the pure N = 2 theory^ The quadratic dependence of the free energy on the dynamically 
generated scale can be viewed as a finite remnant of this quadratic UV divergence. 



2.3 Saddle-point equation 

At large N, the eigenvalue distribution is characterized by a continuous density p(x), defined 
on a interval (-/x, /x). The integral (2.1) is dominated by an eigenvalue distribution determined 
by the saddle-point equation 



fj, ^ 

f dyp(y) I K(x -y)) = 4~ x = ~ 2x ln ( AR ) ■ 

J \x - y J X 



(2.12) 



The eigenvalue density is subject to the normalization condition 



f dy p(y) 



(2.13) 



We now study the solutions to (2.12) in different regimes. 



x This can be seen by setting Ni = 2Ni in eqs. (11), (15) and (18) in [TT]- 
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Small radius 

When AR « 1, one has < A « 1. This is the perturbative regime. In this region, the linear 



we can use the expansion of K(x) in powers of x, 



force represented by 8tt 2 x/X in (2.12) pushes all eigenvalues near x = 0. As a result, fi « 1 and 

(2.14) 



K(x) = -2^(-l) fe C(2A;+l).T 

k=l 



2k+l 



The solution of the eigenvalue density equation (2.12) is then of the form 

/ fen 



\ k=0 I 



(2.15) 



where & max is determined by the perturbative order that one wishes to compute. Then (2.15) 



is a solution to (2.12) with K(x) defi ned b y the sum (2.14) truncated at the same k max . 

Substituting the ansatz (2.15) in (2.12) we find an equation of the form Y,k=o* A k x 2k+1 = 0. 
This gives /c max + 1 conditions for k max + 1 unknowns, a ,ai, ...,a fcmax . Then /z is found by the 
normalization condition. 

Let us explicitly work out the case fc max = 2. Using the formulas, 



t 1 

J dy y 2 " vV-y 2 



\/ti T(n+ \)n 
2(n + 1)! 



l\,,2n+2 



(2.16) 



^ p(y) - = a; 3 ^7rai - -7ra 2y u 2 j + x ^7ra - -na^ 2 - -7ra 2 /i 4 j , (2.17) 

— m 

one gets a linear system of equations for ao, ai, a 2 . Substituting the solution in the normal- 
ization condition leads to a polynomial equation for fi which can be solved in a power series in 
A, 

A 3A 3 C(3) , 25A 4 C(5) 



4vr 2 128vr 6 



2048vr s 



0(A 5 ) 



(2.18) 



The expectation value of the circular Wilson loop (2.5) is then computed from the formula 



W^(C circ i e ) = J dx p{x) 



, 2nx 



(2.19) 



This gives a linear combination of ii(27r/i) and / 2 (27r/i) Bessel functions. Expanding in powers 
of A, we obtain the perturbative series 



\ \2 

W = 1 + - + — — + A 3 
8 192 



/J__5C(3}\ W — 1 7C(3) 35C(5h ^ 

V9216 512vr 4 / V 737280 8192tt 4 8192vr 6 / v ; K J 



This may be compared with explicit evaluation of Feynman diagrams, except that the matrix- 
model calculation is enormously easier and in principle can be pushed to an arbitrary order in 
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A. It is interesting to notice that the expansion can be organized differently, without expanding 
the Bessel functions in A. This would give 



w = - ^ (7r»C(3) - 5C(5)) J 2 (27TA*) + ^ (2 7I a C(3) - 5C(5)) ^(2^/*) 

+ ^^(2^ + (A 4 ) . (2.21) 

Considering the leading term in //, fi « v^A/(27r), the first term in this expansion is the Wilson 
loop VEV in the M = 4 SYM, which is known to resum an infinite series of Feynman diagrams 



without internal vertices [7]. The N = 4 part in (2.20) can be easily recognized as those terms 



that have no £ function coefficients. This part of the series has an infinite radius of convergence. 



It would be very interesting to determine the asympotic behavior of the full series in (2.20). 
It is believed that in gauge theories the number of planar graphs of nth order grows like a 
power of n, as opposed to the n\ total number of Feynman diagrams [?]. This indicates that, if 
the coefficient of the A n term grows as n\, as occurs in the SU{2) cases studied in [6], it would 
be be due to UV renormalons]^] In this case the planar series ( |2.20[ ) would be Borel summable, 
since UV renormalons should not produce poles on the real positive axes of the Borel plane 
(this expectation is consistent with the alternating sign character of the series). 

Large radius 

At large radius the running coupling becomes small and negative. In this regime the eigenvalues 
are strongly repelled from the region around zero. The effective potential for an individual 
eigenvalue looks like the lower curve in fig. [TJ We can imagine two possibilities in this case. If 
the potential wells near the minima are sufficiently deep, the eigenvalue distribution will consist 
of two disconnected patches. Clearly, a two-cut solution cannot be continuously connected to 
the one-cut solution at weak coupling. In this case there will be a large- N phase transition at 
some critical coupling. We found, however, (numerically) that this scenario is not realized in 
M = 2 SYM. The effective potential is always sufficiently shallow, such that the short- distance 
repulsion between the eigenvalues spreads them along a single connected interval for 
any A. For a negative A above some critical value, the eigenvalue density gets a double-hump 
structure with a dip in the middle, shown in fig. [2j 

The extent of the eigenvalue distribution grows rapidly as A -> -oo (exponentially in A 



or linearly in R). The series expansion (2.14) is useless in this regime, because its radius of 



convergence is equal to one. In appendix A we illustrate this point by considering a toy model 



obtained by truncating the series (2.14). The model exhibits a phase transition at a finite R 
and, as a result, the resulting large radius behavior is markedly different from the behavior in 
the actual eigenvalue density in M - 2 theory. 



2 Infrared renormalons arise from the small momentum part of the loop integrals in a certain subset of n-loop 
Feyman diagrams. They should not appear in the present context because the radius of the sphere provides an 
IR cutoff; integrals over momenta are replaced by sums over spin 0, 1/2, and 1 spherical harmonics of S 4 . 
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Figure 2: a) Cases with AR < 1. The eigenvalue distribution expands as A is increased (A = 0.1, 1, 5). 
b) Cases with AR > 1. Eigenvalue density for A = -100, -50, -20. At large AR (i.e. when A approaches 
0~), the eigenvalues expand to a region whose size grows linearly with AR. 



We can solve the saddle-point equation (2.12 ) asymptotically in the limit A ->• , i.e. R -*■ oo, 



where p, ->■ oo. The density will then take a scaling form: p(x) = p(x/p)/p, everywhere apart 
from tiny regions near the endpoints. Using the asymptotic form of the kernel: 



K(x) = 2xln(|a;|e 7 ) + oQJ 



(2.22) 



we can see that in the scaling limit the Hilbert kernel scales away, and we are left with the 
equation 



J~ dy p(y)(x - y) In (|x - y\ e 7 ) = xln AR . 



(2.23) 



Differentiating this equation twice, we then get: 

£ dypjy) 

J x - y 

This equation has a unique normalizable solution: 



. 



(2.24) 



p(x) 



1 



7iy~p 2 - x 2 



(2.25) 



As we will see from the numerics, shown in fig. [3j the analytic solution represents a very 
good approximation to the large radius eigenvalue density, apart from a small region near the 
endpoints, where it violates the boundary conditions. 

Returning to the original equation, one can calculate that 



Jdy 



(x - y) In \x - y\ , ue 

= — = x m — 
,2 2 



K\/p 2 -y 2 



(2.26) 
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P(x) 




5 10 15 20 



Figure 3: Eigenvalue density for A = -15, -10 (solid lines, from top to bottom) and comparison with 
the curves obtained from the analytic formula (|2.25|) (dashed lines). 



Consequently, 

2 

jj = c e ~ x- = c AR , c = 2e ^ 7 . (2.27) 

Let us now compute the expectation value of the circular Wilson loop and the free energy 
in this large radius regime. Using the asymptotic analytic formula (2.25) and (2.19), we obtain 

const / I \ 

W(C chclc ) = 7Tl (27Tfi) » — — - e 2 ^« const ■ exp (47re" 1 " 7 Ai2 - - lnAfl) , AR » 1 . 



(2.28) 

As in the SU(2) case, we obtain a perimeter law. The leading behavior can be interpreted 
as negative mass renormalization Am = -2 e ~ 1_7 A for a heavy probe particle going around the 
circular loop. In terms of the running coupling constant, this is 

In W(C chcle ) « 4vr e -^-^ + ^- + const, A 0". (2.29) 

A 

To compute the free energy, we start with the relation 

dF 
dlnR 

At large R, using the asymptotic density, 

dF 



2N 2 (a 2 ). (2.30) 



dlnR 

Therefore 
at AR » 1. 



-iV 2 /i 2 = -c 2 iV 2 A 2 i? 2 . (2.31) 



F = -— iy 2 A 2 i? 2 , (2.32) 
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Numerical solution 



We start with small values of A, where we know the detailed analytic form of the solution. We 
consider the ansatz 



p( x ) = ( E a k x2k ) \/n 2 -x 2 ■ 



(2.33) 



\k=0 



By taking n sufficiently large, one can achieve any desired accuracy, as long as the solution 
contains only one cut. We have computed the different coefficients and /i numerically for 
n = 1,2,3, ... , and checked that the numerical solution is stable under further increasing n. 



c„ARIfi 




Figure 4: a) [i as a function of AR. b) Comparison between the analytic asymptotic expression (2.27) 
(dashed line) and the numerical value of /u (solid line), where Co = 2e _1 ~ 7 . 

Consider first the interval < A < oo, corresponding to < AR < 1. Fig. [2^, shows the 
eigenvalue density for A = 0.1, 1, 5. The eigenvalue distribution expands as A is increased. The 
same qualitative shape is maintained by further increasing A up to A = oo, where /i has a finite 
value /loo ~ 1.09. As the equation (2.12) has an analytic dependence on the parameter lnAi?, 



-oo. 



there is no discontinuity in crossing from AR < 1 to the region AR > 1. Thus, at A 
the eigenvalue density coincides with A -»• oo. The interesting behavior arises for larger values 
of A. The repulsive force -xlnAR becomes stronger and the eigenvalues are repelled from 
the region of small x. This can be seen in fig. [2]3, which shows the eigenvalue density for 
A = -100, -50, -20. 

The numerical results show that, unlike the toy model of appendix A, the one-cut solution 
exists for all values of A. There is no phase transition. When X/Atc 2 is negative and approaches 
to 0~, in the present case \i begins to increase exponentially with -4tt 2 /A (therefore linearly with 
A_R); the eigenvalue distribution gets stretched out, up to a value ji » 1, with an eigenvalue 
density that, as A -»• 0~, approaches the asymptotic form (2.25). This is shown in fig. [3j which 
displays the eigenvalue density for A = -15, -10 along with the analytic asymptotic eigenvalue 



density (2.25). We see that there is a good agreement between the numerical and analytic 
curves away from the boundary x = ±fi, an agreement which is clearly better for A = -10. 

Fig. shows fi as a function of AR. This figure illustrates in particular that the point 
AR = 1 is not special, as fi(AR) is smooth at this point. It also exhibits the linear law fi ~ AR 
at large R, obtained analytically in the previous subsection. Fig. Hb shows that the slope of the 
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dF 

P = 



dF/dR 




Figure 5: a) dF/dR as function of AR. b) The pressure, p = -dF/dV, as a function of the volume 
V = 87r 2 i2 4 /3. 

curve asymptotically approaches the constant Cq = 2e~ 1 ~ 7 found earlier. Note that the linear 
behavior of /x is a good approximation even in the region AR ~ 1. 

Finally, fig. [5^, shows the derivative of the free energy as a function of AR. At Ai? » 1 it 
exhibits the linear behavior found analytically in the previous subsection. At AR « 1, it has the 
behavior dF/dR ~ N 2 /(2RlnR), as it can be easily deduced from the perturbative formulas. 
We notice that there is a change of sign in d 2 F/dR 2 at the critical radius where d 2 F/dR 2 = 0. 
A natural question is whether this has any thermodynamic implication. The pressure of the 
system is p = -dF/dV, where V is the volume of S 4 , V = 87i 2 R 4 /3. In stable systems, the 
pressure must be positive and, in most common systems, also the compressibility, k = -y^, 
is positive. This is also the case in the present system, as shown in fig. [5]d: the pressure has 
positive sign in the whole interval < AR < oo. The same plot shows that the compressibility 
is always positive. 

In conclusion, we have seen that the one-cut solution exists for all values of AR. We have 
also checked that there are no two-cut solutions at any AR (see also discussion at the end of 
appendix A). 

3 Remarks on superconformal theories 

The simplest Af = 2 superconformal theory (SCYM) is SU(N) SYM with Nf = 2N hypermul- 
tiplets in the fundamental respresentation. The large-A limit of SCYM on S 4 was studied in 
IH [5] , with the idea to get an insight into its putative AdS dual. One may expect that 
in the strong coupling limit it should behave similarly to Af = 4 SYM. Quite surprisingly, the 
strong-coupling limit of the Af = 2 theory is very different from Af = 4. The free energy tends 
to a constant rather than growing logarithmically jl], while the Wilson loop VEV grows as a 
cube of the 't Hooft coupling [3j, in a stark contrast to the exponential growth in the Af = 4 
SYM. 

Below we will comment on the computation of the free energy in Af = 2 SCYM, in Af = 4 SYM 
and in its holographic dual. The holographic computation is a classic result in AdS/CFT \12\ . 
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Yet, we were unable to find a holographic derivation of the In A dependence in the literature, 
and so we will derive this result by reinterpreting the formulas in [12]. We will argue that the 
logarithm arises due to a subtle difference between field-theory and supergravity UV cutoffs. 



3.1 A/" = 4 S YM 

The M = 4 SYM localizes to the Gaussian matrix model, whose exact eigenvalue density is 
given by Wigner semi-circle law: 

p(x) = ^X-4n 2 x 2 . (3.1) 
A 

The free energy then is 

OF 8rr 2 N 2 . 2 . N 2 N 2 

aA = -^H°> = -2A ^ F = -^ lnA ' (3 - 2) 

a result which is easy to obtain directly from the Gaussian matrix integral. 

dF/dX represents the integrated one-point function of the action density (more generally, 
d n F/d\ n is related to the integrated n-point function of the action density). Since dF/dX = 
-Z~ x dZ /dX, the result is independent of the normalization of Z, as long as there is no extra 
A-dependent normalization factor. This requirement is equivalent to saying that the measure 
does not depend on A. 



3.2 Holographic calculation 

At strong coupling, the free energy of M = 4 SYM is identified with the on-shell supergravity 
action on AdS$: 

I = ~ lQ n Q N f d 5 X^(j^R+ —J + /boundary, (3.3) 



r N 

AdS 5 



where Gn is the 5d Newton's constant, L is the radius of the Anti-de-Sitter space, and /boundary 
contains the Gibbons-Hawking surface term and boundary counterterms, necessary to eliminate 
the divergences. 

The metric of AdS 5 , foliated such that the boundary is S 4 , assumes the following form: 



^ 2 = ^ + rW 54 . (3.4) 



The boundary is the large four-sphere at r -> oo. Using the Einstein's equations, R = -20/L 2 , 
the action is calculated to be [T2] : 



4/T r dr r 4 



r drr* 

/ : + counterterms, (3.5) 



3G^L J \/L 2 + r 

where tq is the supergravity cutoff. The counterterms cancel off quartic, quadratic and loga- 
rithmic divergences in this expression. The explicit expressions can be found in [i~2| [T3]. 
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Let us concentrate on the logarithmically divergent term: 



/los = C ln ^ (3 ' 6) 

The counterterm that eliminates this log-divergence is also responsible for the trace anomaly 
in the energy-momentum tensor [T4], so the coefficient in front is directly related to the a = c 
conformal anomaly of M = 4 SYM. Once this coefficient is expressed in terms of N of the 
SU(N) according to the AdS/CFT dictionary: 

it precisely matches with a logarithmic divergence that is also present in the partition function 
of Af = 4 on S A [II]. 

The common lore of the AdS / CFT duality is that one should first eliminate divergences by 
adding counterterms [15] and then compare finite quantities. For the free energy one would 
then obtain a A-independent constant, in a clear contradiction with the field-theory calculation. 
We are going to argue that the reason for this mismatch is a (finite and calculable) difference 
between the field-theory and supergravity counterterms. 

We will relate the supergravity cutoff with the field-theory cutoff by a simple physics argu- 
ment. Consider an SU(N + 1) theory broken down to SU(N) x £7(1) by pulling one brane out 
of the stack that creates AdS and placing it at some radial distance away from the horizon. 
The W-boson, associated with broken symmetry, is described by a string stretching from the 
horizon to the brane. Consider now the W-boson with the largest possible mass. The upper 
bounds on the mass, both in field theory and in supergravity, arise from the UV regularization. 
In field theory, the largest possible mass is obviously A . We can express r through A by 
computing what is the largest possible mass in supergravity. In order to do that we consider 
a heavy probe with the largest attainable mass moving along the big circle of S 4 . On the 
field theory side, the action of the heavy probe is 2-7ri?A . The string action (the area of the 
worldsheet stretching from the horizon to the cutoff surface along the big circle, fig. [6j is 

dr 2nr Ltq 



f f^L = ^ + finite, (3.8) 
from which we conclude that 

27T.RA0 = — - = , 3.9 

a L 

where in the last equality we used the standard AdS/CFT relationship between the radius of 
AdS 5 and the 't Hooft coupling of Af = 4 SYM. 



Combining (3.6), (3.7) and (3.9), we thus express the free energy in terms of the field-theory 
quantities: 

N 2 

^sugra = -— hiA + N 2 \n2irRA . (3.10) 
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The second, infinite piece is renormalized away upon deriving the matrix model by localization)^} 



The finite piece reproduces (3.2). 



3.3 A/"=2SCYM 

Consider now the superconformal M = 2 SU(N) theory with Nf = 2N fundamental massless 
flavors. The partition function is 

C at i -i-r r 9 n i -NY.(^-a 2 +21nH( ai )) 

Z= d^aYlKai-ajfH^ai-a^e 1 >, (3.11) 

i<3 

which leads to the following saddle-point equations: 

f dyp{y) [-±--K(x- y j} = ^x- K(x), (3.12) 



y 

The asymptotic solution at infinite coupling can be found by Fourier transform [3J: 

pM = 2^hf- (3 ' 13) 
The free energy at strong coupling follows then from (a 2 )^ = 1: 

87T 2 iV 2 

F= const + — - — + ... (3.14) 
A 



3 We were informed by V. Pestun that the dependence on R, and consequently on A can be easily recon- 
structed for a generic Af = 2 theory, with the coefficient equal to the a-type anomaly. In particular, in M = 4 



SYM, the Independence is consistent with (3.10). 
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The next correction, as argued in [I], is proportional to In 2 A/ A 2 . We will compute the coefficient 
of proportionality and also the next term in the expansion. 

The asymptotic density extends to infinity, but at any finite A the density is supported on 
a finite interval (-//, /x), where \i is determined by the equation [3]: 



C^JIe^=\, (3.15) 
wherdUC , = 14.60... 

The deviation of the density from its asymptotic form leads to corrections to (a 2 ), which 
are calculated in appendix |Bj 

(«")-!- ^ + ^ + o(f), (3.16) 

where b = 0.4018... The origin of these two terms is easy to understand. The asymptotic density, 
which extends all the way to infinity, underestimates the number of small eigenvalues, the true 
eigenvalue distribution is more compressed. Hence the negative sign of the first correction. 
However, in addition to an overall increase in the density between -/j, and fi, which is the leading 
effect, there is a short-range pile up of the eigenvalues near the endpoints, which produces the 
next term with the positive sign. 



Integrating the free energy with the help of (3.15), which, at this order, amounts to setting 
ji « (2/7r)lnA, we get the final result: 

"8tt 2 



F = const + N 2 



(64tt 2 In 2 A -2^6 In i A + O (In A)) — J + . . . (3.17) 

This is obviously very different from the simple ln(A) behaviour in M = 4, which is exact and is 
valid at any coupling. In the SCYM, the free energy is a non-trivial function of A, and behaves 
logarithmically only in the weak-coupling limit. 

4 Remarks on Af = 2* theory 

This theory is obtained by deformation of M - 4 SYM by adding a mass^ M term for the 
adjoint hypermultiplet. In this case the equation for the eigenvalues is 

/ dyp(y) (—-K(x-y) + \{K(x -y + M)+K(x-y-M))) = ^ x . (4.1) 
J \x-y 2 / A 

Now the eigenvalue density depends on a two-parameter space Xq,M. Some corners of this 
two-parameter space were studied in jB]. For small M, it approaches the eigenvalue density of 
Af = 4 SYM plus 0(M 2 ) corrections, 



p(x).v^^^(l-M 2 ), A i= }^-= , (4-2) 
Ao 2tiV1 - M 2 



An analytic expression for the constant C in terms of an infinite sum is given in [3]. 

5 Here the radius of the sphere is absorbed into the definition of the mass of the hypermultiplet: M -* MR. 
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for M « 1 and Ao » M 2 . In particular, one finds W(C chcle ) w exp (1 + M 2 /2)] and 

N 2 

F w — — (1 + M 2 ) In A , M«l, A »M 2 . (4.3) 
For large M and fixed M e A » = A, the theory flows to the J\f = 2 SYM theory studied 

47T 2 

in section 2, with renormalized coupling A defined by AR = e a . To show this explicitly, we 
take the large M limit, where one can approximate 

2K(x -y)- K(x -y + M)- K(x - y - M) = 2K(x -y)- 2K'(M) (x - y) + 0(1/ M 2 ) (4.4) 

Thus, neglecting 0(1/M 2 ) contributions, we have 

/ dyp(y) (—-K(x-y) + K'{M)(x -y)) = ^x. (4.5) 
J \x-y } A 

-i" 

Using 

J dyp(y)K'(M)(x - y) = xK'(M) , (4.6) 
-n 



and K'(M) » 21nMe 1+ ^ + 0(1/M 2 ), we are left with the saddle-point equation Q2.12D for the 
Af = 2 pure SYM theory. 



Returning to the M = 2* theory, by solving the saddle-point equation (4.1) numerically 
one can also explore other regions of the two-parameter space Ao, M, away from the corners 
represented by the M = 4 and M = 2 SYM theories. One novel feature with respect to the M = 2 
theory studied in section 2 is that, for a wide range of parameters, eigenvalues accumulate near 
x = and near x = ±/i at the same time. 

This is illustrated by fig. [7j It exhibits a hill at x = and two more at x = ±/x. As M 
is further increased, the hill at x = diminishes, and \x increases, approaching the eigenvalue 
density of J\f = 2 theory studied in the previous section. 

We have also checked that in J\f = 2* SYM theory the one-cut solution exists in all regions 
of the parameter space and that there are no two-cut solutions in any region of the parameter 
space. 



5 Instanton contributions 

Instanton contributions are usually believed to be negligible at large N, because of the expo- 
nential factor exp(-87r 2 N / X) in the instanton weight. It is possible, however, that the instanton 
moduli integration overcomes this exponential suppression making the instanton contribution 
formally divergent in the large- limit. This phenomenon leads to an instanton- induced phase 
transition [9]. It was shown in [3] that such a transition does not happen in the superconfor- 
mal M = 2 theory discussed in sec. |3.3| and that the instantons always remain exponentially 
suppressed in M = 2 SCYM. Here we analyze the instanton weight in the pure J\f = 2 SYM. 
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Figure 7: Eigenvalue density for A = 107T 2 , M = 9 in Af = 2* SYM theory. 



The expression for the instanton factor in pure M = 2 SYM theory can be obtained by 
taking a suitable limit of the instanton weight in M = 2* theory. When regarding the M = 2* 
theory as a UV regularization of the pure M = 2 SYM some care is necessary in picking the 
same renormalization scheme in all parts of the localization partition functions. By that we 
mean the use of the same definition for the running coupling, which is a priori ambiguous due 
to the choice of the subtraction point. As in sec. |4| we adopt the following definition: 



4tt 2 



-lnM-1-7. 



(5.1) 



where Ao = g 2 N represents the bare coupling of M = 2* and as before A denotes the running 
coupling constant of the M = 2 theory. We emphasize that once we have used this renormal- 
ization prescription in evaluating the perturbative part of the partition function, the definition 
of the running coupling in the instanton weight is not a matter of choice any more. We should 
use the same definition. In particular, the dynamically generated scale, 



4vr 2 



-InAR, 



is then given by 



AR = Me 



-+1+7 



(5.2) 



(5.3) 



The pure M = 2 SYM without matter multiplets arises in the limit M 00, Ao -> with 
A fixed. The one-instanton contribution to the partition function of the M = 2* theory can be 
found e.g. from the expression given in [16] (which is a particular case of the Nekrasov partition 
function [17]): 



Z} nst = - c 



2 N 

i=\ j*i 



(ai-dj+i) - M 2 
(oj - a,j) {a t - a,j + 2i) 



(5.4) 
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Here 8 Q is the theta-angle of the J\f = 2* theory. Another possible representation is through a 
contour integral: 



Z** - 2M2 T - fl (Z "^ )2 - M2 (5 5) 



where the contour of integration encirlces the poles at a,- + i counterclockwise. 



Upon taking the limit M -> oo, \ -*■ with A fixed, the instanton weight (5.5) becomes 



y^N+ie r dz SL 



= e -(^—j— / _ n - — Vt> (5 - 6) 

,=i (2 - + 1 



with the renormalized coupling given by (5.1). Remarkably, for odd N, the theta-angle also 
gets renormalized: 

9 = 6 + irN. (5.7) 
We recall that the running coupling A becomes negative at AR > 1. In such a situation, the 

8 2 

factor e ~ does not suppress the instanton contribution, but, on the contrary, enhances it! 
One might naively conclude that instantons become very important at large N if AR > 1. 

To clarify this important point, we need to compute Z{ nst . In the large- N limit, the instanton 
weight becomes 

Z'r = A^ +2+2 ") N+ld r^ e -"*-«, (5.8) 

J 7T 

where the z integration can now be taken along the real axis, and 



+1-1 



•S'instC 2 ) = J dx p(x) In ((2 - x) 2 + l) . (5.9) 

-/< 

The z- integral is of the saddle-point type, with the single saddle-point at z = 0, so we find that 

^inst „ Q -NS cf! ( 510 ) 

with 

S cS = -21nA J R + 2 + 2 7 + ^ rfxp(x) In (x 2 + l) . (5.11) 

-/< 

The instantons are suppressed if and only if S e s > for any value of AR, which is not at all 
obvious as the first term becomes negative for AR > 1. We have checked that the full action 
is nevertheless positive for any AR, by computing the action analytically in the two limiting 
cases of AR -*■ and AR -*■ 00 and numerically at intermediate Ai?'s. 

At small AR, the integral in (5.11) can be estimated as (ln(l + a; 2 )) ~ (x 2 ) ~ ~ A ~ 
-1/lnAi?. The action thus is large and positive. Obviously, as AR grows, the action will 
diminish (at first logarithmically). The question is whether it crosses zero or not. To answer 
this question we can study the opposite limit of very large AR. 

8 ^JV 

When AR ->■ 00, the bare instanton weight, A -> 0~ and e ~ ->■ 00 , which would naively 
imply that the one-instanton contribution blows up at N -> 00, but we also need to take into 
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account the last term, the matrix- model average of ln(x 2 + 1). Since we know the solution of 
the matrix model analytically at AR -*■ oo, we can take this solution and calculate the integral. 



Inserting p(x) from (2.25), and taking into account that \x » 1, we get: 

J dxp(x) ln(l + x 2 ) = 2 In i-\J p? + 1 + - j 



21n^ + - + 



(5.12) 



Using the explicit expression for \i from (2.27), we finally get 



, 1+7 



cff 



AR 



+ Ol 



(ARy 



(5.13) 



The instanton action becomes very small at large AR, but remains positive. 

We have checked numerically that the instanton action is a monotonic function of AR and 
consequently never becomes smaller than zero. In conclusion, in the large N limit, the one- 
instanton contribution is negligible. 



6 Conclusions 

We found that Wilson loops in M = 2 SYM on a large four-sphere obey a perimeter law with a 
precisely calculable coefficient. The free energy scales as R? . It would be really interesting to 
apply localization techniques to confining theories, perhaps by breaking M = 2 supersymmetry 
to jV = 1 as in [8]. The supersymmetry breaking unfortunately ruins localization, because 
the superpotential term that one should add to the Lagrangian does not commute with the 
supercharge used in the localization procedure. 

The large- N superconformal Af = 2 theory with Nf = 2N fundamental hypermultiplets turns 
out to be very different from Af = 4 SYM at strong coupling [2J |3l HJ |5] . The dual string theory, 
if such exists, should thus remain strongly coupled, even when the 't Hooft coupling is large, 
in contradistinction to the Af = 4 case, where at large 't Hooft coupling the dual string theory 
behaves semiclassically at large A. 
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A Toy model with phase transition 

It is instructive to consider a toy model, obtained by truncating the expansion of K(x) to the 
first term. We replace 

K(x) -+ 2((3)x 3 . (A.l) 



The equation for the eigenvalue density (2.12) can then be solved exactly by the following 
ansatz: 

p(x) = (n + cx 2 )V '/i 2 - x 2 , (A.2) 
where no, c, \i are parameters to be determined. The normalization condition gives 



2 cfi 2 x 
no = — - -j- ■ A.3 



In order to determine w and c we use the formulas 

M 1 7T 

y Gfo/(n + °y 2 )^/ V 2 ~ y 2 ( x _ y ^ = 7in ^ x + ^ cx{2x 2 - « 2 ) , 

/I 7T 7T 

dy{riQ + cy 2 )\f jj? - y 2 {x - y) z = — noXfx 2 (Ax 2 + 3/i 2 ) + — cx/i 4 (2x 2 + 3/x 2 ) 
8 16 



Then the saddle-point equation (2.12) becomes 



n + 1 c(2x 2 - « 2 ) - C(3)| w 2 (n (4x 2 + 3/i 2 ) + | /i 2 (2x 2 + 3 W 2 )) = y . (A.4) 



This leads to 

» 2 C(3) c= 2C(3) 

7T/i 2 27T ' 7T 

and \x 2 = z is the positive real root of the quartic equation 



no = ~^-, c-i^-, (A.5) 



A 2 16 



(A.6) 



Using these expressions one can compute the eigenvalue density for different values of A. 
For small A > 0, the first few terms in the perturbative expansion are as follows 

+ ™ + (A.7) 



4tt 2 128vr 6 16384tt 



This shows that, in this weak coupling regime, all the eigenvalues x are confi ned in a small 
region 0(\/A), so that the term K(x) = 0(x 3 ) is small. Comparing with ( 2.18| , we see that in 
this regime the difference between the toy model and the actual J\f = 2 theory begins at order 
0(A 4 ). 
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When A ->• ±00, \i approaches the asymptotic value, 











A=±oo \ 


(>/3 1 ) C(3) 



0.72 



and 



n 



A=±oo 



2 

7T 



3-1 




+ 1 C(3) *0.75 . 



(A.i 



(A.9) 



As A is increased from -00, the eigenvalues begin to accumulate near some x + 0. The 
eigenvalue density at x = begins to decrease until it vanishes at the critical value where 
n = 0. This occurs at 



A, 



71" 



A/03) 



>/C(3) ' 



(A.10) 



For A cr < A < 0, the system no longer admits the one-cut solution (|A.2|). Now the eigenvalue 

(A.ll) 



density has support in two disconnected regions. The new solution is 

p(x) = n\\J (a 2 - x 2 )(x 2 - b 2 )x 2 . 



To perform the integrals in (2.12) we now use 



f dy^(a 2 -y 2 )(y 2 -b 2 )y 2 = ^(a 2 -6 2 ) 2 



x(2x 2 -a 2 -b 2 ) 



f dyyj (a 2 - y 2 ){y 2 - b 2 )y 2 -^— = 
J x - y 

f dyJ{a 2 - y 2 ){y 2 - b 2 )y 2 (x - y) 3 = £ x(a 2 - b 2 ) 2 (2x 2 + 3a 2 + 3b 2 ) (A.12) 
J 16 



Solving the normalization condition and (2.12) we find 

r2 



a = — 



+ 



1 



b 2 = - 



1 



2C(3) 



(A. 13) 



A C(3) ' VC(3) ' " A C(3) > /c(3) ' 

Since a 2 and b 2 must be positive, the solution exists only for A cr < A < 0. Thus this solution 
appears precisely at the point the one-cut solution (A. 2) ceases to exist. 

Let us now compute the free energy before and after the phase transition. The free energy 
is given by the formula 



l7T 2 N 2 

A 



J dxp(x)x 2 -^N 2 J dx J dyp(x)p(y)(\n(x-y) 2 + \nH 2 (x -y)) (A. 14) 



In the toy model, 

\nH 2 (x) ^-((3)x 4 . (A.15) 

To understand the character of the phase transition, we only need the derivative with respect 
to A. On the solution, we have 



dF _ OF 

~d\~~d\ 



57T 



'■N 2 



A 2 
21 



J dxp(x) 



x 



(A.16) 



In the phase I lying in the regions < A < oo and -oo < A < A cr , we have 

f ^ Pl (a;)x 2 = ^(4 + C(3)/i 4 ) . (A.17) 



where /i = /i(A) is defined by the positive real root of (A. 6). In phase II, lying in the region 
A cr < A < 0, we have 

J dxpn(x)x 2 = -J^ ■ ( A - 18 ) 

Therefore 

dFi 8n 2 N 2 fi\ A A/oX 4X dF n 8n 2 N 2 tt 2 _ 

By further differentiation, we can compute second and third derivatives of the free energy. At 
the critical point A cr 



7T 


we find 






7c(3) 








dFi 


dF u 


8iVVC(3) 




~dX ' 


dX 






d 2 Fj 


d 2 F u 


24N 2 ((3) 




d\ 2 


dX 2 


7T 4 




d?F l 


98N 2 ((3)I d 3 F n 


96A^ 2 C(3)i 


dX 3 




7T 6 ' dX 3 


7T 6 



(A.20) 

Thus there is a discontinuity in the third derivative of the free energy: the phase transition is 
of the third kind. 

Next, we compute the expectation value of the circular Wilson loop in both phases I and 
II. We use 

W(C ciTcle ) = J dxp(x)e 2 ™ . (A.21) 
Hence, in phase I, with eigenvalue density (A.2), we find 

^i(C circlc ) = (7r/a(27r^) (3//C(3) + 4) - 6^ 3 C(3)/ 2 (27r/i)) . (A.22) 

At small A, we find the expansion 

^i(C circlc ) = l + - + — + X 3 l— - + O (A 4 ) . (A.23) 

J 8 192 V9216 512tt 4 / y J K J 

Up to this order, there is no difference between the toy model and J\f = 2 SYM. The difference 
starts at the A 4 terms (cf. (2.20)). As A -> oo, Wi(C c \ TC \ e ) approaches a finite constant, since \i 
approaches the asymptotic value (A. 8). 

In phase II, we have 

Wn(C circle ) = ^® f dx^{a 2 -x 2 ){x 2 -b 2 )x 2 e 2 ™ . (A. 24) 

7T J 
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We consider two limits: 

a) When A ~ A cr , then b ~ 0, the integral can be computed with the result 

a 2 fCV\ 2 
W n (C circle ) = -^/(27fal 1 (2an) - 3I 2 (2an)) , a 2 = -j== . (A.25) 



This of course agrees with (A. 22) at the critical point since at this point the eigenvalue densities 
pi and pn are the same. 

b) When A -»• we have that a,b -> oo. By rescaling the integration variable, the integral 
then collapses to a point and can be easily computed. We find 

W u (C circle ) « - exp[27r 2 /V-X(3)] • (A.26) 



7T 



In terms of AR, this is 



W n (C ckcle )«- eVS«. (A.27) 

7T 

Thus the behavior of the Wilson loop in this toy model completely deviates from the behavior 
of the Wilson loop in the original SYM system of section 2, where H^(C circ ie) obeyed a perimeter 
law. 

One may also consider a more general toy model obtained by replacing K -»■ K e in (2.12) 
where 

K e (x) = e~ 2 x (ip(l + iex) + ip(l - iex) + 27) (A.28) 

This interpolates between the toy model ( |A.l ), recovered by taking e -> - since K e (x) = 
2((3)x 3 + 0(e 2 ) and the actual M = 2 matrix model, corresponding to e -*■ 1. For sufficiently 
small e the model still undergoes a phase transition towards a phase represented by the two-cut 
solution. A natural question is for which value of e the two-cut solution disappears. One obtains 
the following picture. Numerically, for sufficiently small e, one finds two-cut solutions for A in 
an interval (A cr ,0~), with a similar A cr . However, as e is increased, one reaches a critical value 
of e ~ 0.3 beyond which the two-cut solution does not exist. In particular, it does not exist in 
the case of interest, e = 1. 



B Free energy in Af = 2 SCYM 



Here we present the derivation of eq. (3.16). The starting point is an approximate solution for 
the density obtained in [3]: 



p{x) = J 



with 



1 2 sinh 2 



2~i 



F(u) 



e P(w) 



(B.1) 



coshw cosh a; 

r n e -i^ n F(-iv n ) {y n - u ) (ui - iv Q ) 



n=l 



(u n + uq) (u + iv Q ) (u + iv n ) 



(B.2) 
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where 

^ * = Xuo A ' ^ ' 

The function 



G(w) = 7 . V W (B.4) 

is analytic in the lower half-plane, and has simple poles at uj = %v ny 

v n = 7r(n+-\, (B.5) 

with residues 

r n = res G(u) = -— p— — . (B.6) 

w =^n ^/n [n + ij r (n + 1) 

The constant a has a numerical value a = 1.0232... Its representation through a combination of 
infinite sums can be found in [3]. 

The approximate density was obtained with the help of the Wiener-Hopf method and is 
pretty accurate on the whole interval (-/i, /i), except for a small vicinity of the left endpoint, 
-fi. In particular, it violates the boundary condition p(-yu) = and thus should be used with 
caution. A safe strategy is to use it only for positive x. 



The first term in (B.2 ) is the Fourier transform of the asymptotic density (3.13 ). The second 
term is the first correction to the density due to the finiteness of /i. This term is small when 
fi is large, but affects the density evenly on the whole interval (-/i, /i). The last term, on 
the contrary, is rapidly oscillating, and consequently is only important in a small vicinity of 
the endpoints. The purpose of this term is to subtract the poles in the lower half-plane, thus 
ensuring that p(x) = for x > fi. 

When using the approximate solution above, all averages should be computed as integrals 
over positive x. In particular, 



For the first two terms in (B.2), which are manifestly even functions of u, we can use the 
Sokhotski-Plemelj formula, giving in this case the expected p"(0). This trick does not work 
with the last term, because the principal value integral does not vanish. Instead, the last term 
should be evaluated by closing the contour in the upper half-plane and computing the residues: 

{a ) = l-F(0) + V2/i^ , (B.8) 

n=l V n \V n + Vq) 

where we have only kept the leading terms at \i -*■ oo. Taking into account that at large \x 
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we get the formula in the text, eq. (3.16), with 



2jr g r n n 
n n=i (2n + 1)2 (n + 1) 



a(2n + l) 
8^ 



0.4018... (B.10) 
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